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Abstract. The bispectrum of primordial curvature perturbations in the squeezed configu- 
ration, in which one wavenumber, ks, is much smaller than the other two, fcs ^ fci ~ k2, 
plays a special role in constraining the physics of inflation. In this paper we study a new 
phenomenological signature in the squeezed-limit bispectrum: namely, the amplitude of the 
squeezed-limit bispectrum depends on an angle between ki and such that B(^{ki, k2,k^) — )• 
2'^j^CLPL{k.i • k3)P^(A:i)P|j(/c3), where Pl are the Legendre polynomials. While cq is related 
to the usual local-form /nl parameter as cq = 6/nl/5, the higher-multipole coefficients, ci, 
C2, etc., have not been constrained by the data. Primordial curvature perturbations sourced 
by large-scale magnetic fields generate non- vanishing cq, ci, and C2. Inflation models whose 
action contains a term like I{(j)f'F'^ generate C2 = co/2. A recently proposed "solid inflation" 
model generates C2 S> cq. A cosmic- variance-limited experiment measuring temperature 
anisotropy of the cosmic microwave background up to ^max = 2000 is able to measure these 
coefficients down to 5cq = 4.4, 8ci = 61, and 5c2 = 13 (68% CL). We also find that cq and 
ci, and Co and C2, are nearly uncorrelated. Measurements of these coefficients will open up a 
new window into the physics of inflation such as the existence of vector fields during inflation 
or non-trivial symmetry structure of inflaton fields. Finally, we show that the original form 
of the Suyama-Yamaguchi inequality does not apply to the case involving higher-spin fields, 
but a generalized form does. 
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1 Introduction 

Convincing detection of the so-called "local-form" three-point correlation function (bispec- 
trum) of primordial curvature perturbations from inflation would rule out all single-field 
inflation models [1, 2], provided that an initial quantum state of the curvature perturbation 
is in a preferred state called the Bunch-Davies state [3, 4] and that the curvature perturbation 
does not evolve outside the horizon due to a non-attr actor solution [5, 6].^ 

The curvature perturbation, C, is defined as a trace part of space-space components of 
the metric perturbation, 6gij = a?{t)e'^^5ij, in a uniform density gauge. The bispectrum of Q 
is defined as (CkiCk2Ck3) = (27r)^(5^^) (ki + k2 + ka) B(^{ki, k2, k^), and the local-form bispec- 
trum is defined as Bi^{kiM.h) = |/nl [Pc{ki)P^{k2) + Pc{k2)Pc{hi) + Pi^{k^)P^{ki)] (e.g., 
[7]), where P(;{k) oc k^''~'^ is the power spectrum of the primordial curvature perturbation 
with Us = 0.96 lb 0.01 [8-10]. This means that the local-form bispectrum has the largest 
amplitude in the so-called squeezed configuration, in which the smallest wavenumber, k^, is 
much smaller than the other two, i.e., fcs ^ /ci ~ k2 [11]. In this limit the local-form bis- 
pectrum is given by — )• ^/nl-Pc(^i)-^c(^3)' ^^"^ attractor single-field inflation models 
with a Bunch-Davies initial state give /nl = ^(1 — ?^s) [1, 12]. 

The current best limit on /nl is /nl = 37 it 20 (68% CL), which was obtained from the 
Wilkinson Microwave Anisotropy Probe ( WMAP) 9-year data with the expected IS W-lensing 
bias removed [13]. The forthcoming Planck data are expected to reduce the error bar by a 
factor of four [7]. 

^Also see workshop summaries of "Critical Tests of Inflation Using Non-Gaussianity" in http://www.mpa- 
ga rching.mpg.de/~komatsu /meetings/ng2012/. 
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Figure 1. Absolute values of the shape hmction of L = 1, (/ci/c2fc3)^5'i (left panel), and that of L = 2, 
(fcifc2fc3)^5'2 (right panel). We restrict the plot range to < k2 < ki and \ki — k2\ < k^ < ki + k2 
for symmetry and the triangular condition. The shape of L = 2 peaks at the squeezed configuration, 
ks/ki <C 1 and fc2/fci ~ 1, whereas the shape of L = 1 is suppressed at the squeezed configuration. 



If the Planck collaboration finds evidence for /nl, or the lack thereof, what is next? 
Measuring the local-form four-point function (trispectrum) [14, 15] to check the so-called 
Suyama-Yamaguchi inequality between the amplitude of the local- form trispectrum and /nl, 
i.e., tnl > (6/nl/5)^ [16-21], would be an important next step to understand the nature of 
sources of non-Gaussianity (or the absence thereof). We shall discuss the Suyama-Yamaguchi 
inequality within the context of higher-spin fields in Sec. 4. 

Can we learn more about sources of non-Gaussianity by further scrutinizing the behavior 
of the bispectrum in the squeezed configuration? The answer is yes, and this is the main goal 
of this paper. Namely, in this paper, we shall investigate phenomenological consequences of 
the following new parametrization of the bispectrum of primordial curvature perturbations: 

B^ih,k2,h) = Y,clPlC^i ■U)Pdh)Pdk2) + {2 perm) , (1.1) 

L 

where Pl{ij) is the usual Legendre polynomials, i.e., i-b(/u) = 1, -Pi(a*) = A*, and -P2(/^) = 
^(3//^ — 1). Here, cq is equal to 6/nl/5.^ 

Why consider cl with L > 1? These coefficients appear to be sensitive to the existence 
of vector fields. For example, curvature perturbations sourced by primordial magnetic fields 
produce non-zero ci and C2 [22, 23]. Curvature perturbations sourced by a jI{4>)'^F'^ term in 
Lagrangian produce C2 = cq/2 [24, 25]. These coefficients are also sensitive to the existence of 
a non-trivial realization of S0(3) rotational symmetry during inflation: a recently proposed 
"solid inflation" model produces C2 ^ cq [26]. While the second-order effects in General 
Relativity also induce non-trivial angular dependence in the bispectrum, it disappears in the 
squeezed limit [27], and thus will not be considered in this paper. 

While we assume that the coefficients, c^,, do not depend on wavenumbers, it is entirely 
possible that they do. There are various ways in which cq = 6/nl/5 depends on wavenumbers 

^Note that, due to symmetry, the ci term as well as any odd L terms vanish in the exact squeezed limit, 
i.e., limfeg^o B(;(ki, fc2, fcs) = 2[co + C2P2(mi3) + C4P4(/ii3) + . . .]P(;{ki)P(;(k3), where ^13 = ki -ka. As a result, 
when we analyze the CMB data, the error bar of Ci is much bigger than the error bars of co and C2, as we 
shall show in Sec 3.3. On the other hand, the error bars of cq and C2 are expected to be comparable: we shall 
show that they are related by Sc2 ~ 3(5co in Sec 3.3. 
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[3, 4, 28-32]. Particularly interesting possibilities are strongly infrared-divergent ci and C2, 
which can naturally give rise to dipoler (i.e., hemispherical) and quadrupolar modulations, 
respectively, of the observed power spectrum in our sky [33] . 

In Sec. 2, we shall briefly review these three scenarios to motivate our choice of parametriza- 
tion given in Eq. (1.1). Specifically, we review non-Gaussianities generated from: (1) large- 
scale magnetic fields after inflation in Sec 2.1; (2) a vector field coupled to the inflaton field, 
(/>, through a dilaton-like coupling P{(j))F'^ in Sec 2.2; and (3) solid inflation, in which the 
inflaton field is a part of a vector multiplet, { c/)^ , (/)^ , (/)^ } , in Sec 2.3. 

What do the shapes of L = 1 and L = 2 terms look like? Using the triangular condition 
of three wavevectors, ki -|- k2 -|- ks = 0, the cosine between wavenumbers can be written 
as, e.g., kl • k2 = (A:| — kf — /c|)/(2A;iA;2). Then, for a scale-invariant power spectrum of C) 
P({k) = ^As, Eq. (1.1) can be re-written as 

B^ikiMM) = {2Ti^AsfY.''LSL{kiMM) , (1-2) 

L 

where 



5'o(A;i,/c2, /cs) = T373 + 2 perm ) , (1.3) 
Si{ki, k2, k^) = [i^tu + 2 perm J - ( — ^ + 5 perm ) , (1.4) 

/ 3/(;4 \ / 3^2 

S2{ki,k2, fcg) = (^^^ + 2 perm J - (^^^ + 5 perm 

+ (^77l7- + 5perm) + ( — i-^+2perm) . (1.5) 



Sklk2 J \ikfkl 

In Figure 1, we show the shape functions Si (left panel) and S2 (right one). The right 
panel shows that the bispectrum for L = 2 peaks in the squeezed limit (fes <^ ki ^ k2), 
whereas the left panel shows that the bispectrum for L = 1 is suppressed in the squeezed 
limit (also see footnote 2). 

The rest of this paper is organized as follows. In Sec. 3, we derive both the full-sky and 
flat-sky formulae of the bispectrum of temperature anisotropies of the cosmic microwave 
background (CMB) induced by the angle-dependent bispectrum given in Eq. (1.1), and 
analyze their behaviors. We also estimate the error bars of cl for L = 0, 1, and 2, ex- 
pected for a cosmic-variance-limited CMB experiment measuring temperature anisotropy up 
to 

^max — 2000. In Sec. 4, we revisit the Suyama-Yamaguchi inequality within the context 
of higher-spin flelds such as those discussed in this paper. We conclude in Sec. 5. In Ap- 
pendix A, we discuss the precision of the fiat-sky approximation. In Appendix B, we derive 
the CMB bispectrum in the Sachs- Wolfe limit (in which the temperature anisotropy is given 
by 6T/T = — (^/5). In Appendix C, we present the full Fisher matrix for cq, ci, and C2- 

Throughout this paper, we adopt the following convention for Fourier transformation 
of an arbitrary function, /(x): /(x) = J ^|^/(k)e*'^'^. 

2 Theoretical motivation 

2.1 Helical and non-helical magnetic fields 

Astrophysical observations suggest the existence of magnetic fields on the order of 10~^ G in 
galaxies and cluster of galaxies [34-37]. There is also indirect evidence for the existence of 
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magnetic fields on the order of 10^^^ — 10~^^ G in the inter-galactic medium (IGM) [38-41]. 

There is yet no compelling model for how these vector fields can be generated during 
inflation, as the existing models suffer from strong backreaction or strong coupling problems 
[24, 42-46] (see the more detailed discussion in the next subsection). Here we simply assume 
that a magnetic field has been generated,^ and study its impact on the primordial pertur- 
bations during the radiation era and recombination. In the next subsection, we discuss the 
additional signatures that take place if vector fields are coupled to the inflaton field. 

A large amount of literature exist in the studies of effects of vector fields on CMB 
anisotropies and the large-scale structure of the universe. See, e.g., [47-60] for effects on 
the two-point correlation functions, and [22, 23, 61-67] for those on higher-order correlation 
functions. 

Let us assume that super-horizon vector perturbations were produced during infiation, 
and they generated large-scale magnetic fields. The anisotropic stress of this magnetic field 
sources the growth of curvature perturbation via Einstein's field equations during the radia- 
tion era. However, after the decoupling of neutrinos at a few MeV, the magnetic anisotropic 
stress is compensated by the neutrino anisotropic stress, and the curvature perturbation on 
super-horizon scales becomes a constant. This constant curvature perturbation survives till 
the recombination epoch, and seeds additional CMB anisotropies. The solution of curva- 
ture perturbations on super-horizon scales is determined by the traceless projection of the 
magnetic anisotropic stress as [53] 

&.0.91n(^) (i4.-i.'.)^/|i;B-m(k-k'), (2.1) 

where tb and denote the conformal time of the generation of magnetic fields and that 
of the decoupling of neutrinos, respectively, and p-yfi is the present-day value of the photon 
energy density. This equation shows that, even under the assumption that the magnetic field 
itself is a Gaussian variable, the curvature perturbations become highly non-Gaussian. 

Angular dependence of the power spectrum and bispectrum arises due to the spin-1 
nature of magnetic fields. The magnetic field vector is transverse, k^Bi = 0, and thus it is 
expanded using the spin-1 polarization vector, ej'^^(k): i?i(x) = (27r)~^ J dPh Bi(k)e^^'^ = 

(27r)~^ f d^k -B'^(k)e|'^^ (k)e*'*''', where a denotes two circular polarization states. 

Then, the power spectrum of magnetic fields can be decomposed into the "non-helical," 
Psik), and "helical," Peik), components as [68] 

{B,{k)B,{k')) = i^<^(3)(k + k') [PB{k)-aP,s{k)]e'f\k)e^-''\k) 



(2vr) 



-(^(3)(k + kO [{6ij - kik)PB{k) + ie,jiUPB{k)\ , (2.2) 



^While we use the term "magnetic field" and "electric field" here and in the next subsection, these fields 
are not necessarily the usual electromagnetic fields. For the discussion in this subsection, it is sufBcient to 
have some vector field whose anisotropic stress decays as Tj — ^SjTj^ oc on super-horizon scales. On the 
other hand, the anisotropic stress on super-horizon scales is constant during inflation (disregarding slow-roll 
corrections) for the case we discuss in the next subsection. The important feature of these models is that the 
anisotropic stress scales with a in the same way as the isotropic pressure dominating the universe: for the 
former case, it scales in the same way the radiation pressure does, and for the latter it scales in the same way 
the inflaton pressure does. 
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where Ciji is the antisymmetric tensor normahzed as £123 = 1. These power spectra are 
defined as 



(i?+(k)i?+(k')) + (i?-(k)i?-(k')) = (27r)3Ps(fc)5(3)(k + k'), (2.3) 
-(S+(k)B+(k')) + (S-(k)S-(k')) = {2TrfPB{k)6^^Hk + \^'). (2.4) 

Note that the overall sign of the definition of these spectra depends on the choice of the 
polarization vector. Here, we adopt the definition of Ref. [23]. 

Using these power spectra, one can write the angle dependence of the bispectrum of 
curvature perturbations as [23] 

/I 2 

B(;(ki,k2,k3) OC PB{h)PB{kl)PB{k2) ( 2/"l2 + ^23 + ^31 " 3 " ^*12)"23/ti31 
- PB{k*)PB{kl)PB{k2) ( ^23/J'31 " ^Ail2 



3' 

+ (1^3, 2^1, 3 ^ 2) + (1 ^ 2, 2^3, 3^1), (2.5) 

where fiab = k^ • k^ and A;* denotes some pivot wavenumber.^ 

Eq. (2.5) clearly shows that helical and non-helical magnetic fields generate L = 1 and 
L = 2 angular dependence in the bispectrum of curvature perturbations. If the magnetic 
field was generated at a GUT scale, i.e., t^/tb ~ 10^'^, with nearly scale-invariant spectra of 
Pb and Pq, the Legendre coefficients in Eq. (1.1) are related to the amplitudes of non- helical 
and helical magnetic fields smoothed on 1 Mpc as [23] 

— 10- (i^)'. . . -0. (^)^ {^)\ . . 14... (2.0) 

Therefore, if inflation creates i?iMpc ~ 3nG and ^iMpc ~ InG, which are consistent with 
the current observational limits, we may have negative and non- vanishing ci and C2; namely, 
Cl ~ —8 and C2 ~ —2. 

2.2 P{(I))F^ model 

Vector flelds with the standard Maxwell — kinetic term are not produced by the ex- 
pansion of the universe and, if generated by some other source, they are rapidly diluted 
away. This poses a challenge to models of primordial magnetogenesis and of vector fields 
during inflation. Vector fields during inflation can result in broken statistical isotropy of the 
primordial perturbations, which will be probed by the forthcoming Planck data [70, 71]. 
Vector fields with a kinetic term given by 

£ = -^F^ (2.7) 

can instead be produced during infiation if / (t) has an appropriate time dependence [72] . It 
is convenient to define the "electric" and "magnetic" components 

Ei = , Bi = -^€ijkdjAk , (2.8) 



''Note that parity-odd terms, which are proportional to PgPB and Pg, do not appear in Eq. (2.5), as is 
a scalar. On the other hand, the bispectrum involving vector or tensor perturbations may contain parity-odd 
terms, which yield the CMB temperature auto-bispectrum with l\ -\- 12 + 13 = odd [23, 69]. 
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where the primes denote derivatives with respect to the conformal time, and a is the scale 
factor of the universe. In terms of these components, the physical energy density in the 
vector field assumes the conventional expression, pA = 2 — ^• 

If / oc a", with n = +2 or n = — 3, the magnetic modes are generated during inflation 
with a scale invariant and frozen spectrum outside the horizon. Rather than assuming that / 
is an external function, one can obtain the required time dependence by assuming that I is a 
function of the inflaton (j), with a functional form related to the inflaton potential by [72, 73] 



I = Iq exp 



n d(j) 



{I) (X a", (2.9) 



where e is the usual slow-roll parameter, e = (f^) ' '^i^h Mp = l/\/87rG denoting the 
reduced Planck mass. 

Some recent work studied whether this coupling can result in visible cross-correlations 
between primordial perturbations and large-scale magnetic fields [58-60, 74-77]. This is not 
trivial to realize, as the n = —3 choice results in too large an energy density in the electric 
modes [42, 43], while n = +2 leads to too large an electromagnetic coupling constant during 
inflation [24, 43]. ^ 

For these reasons, we prefer not to identify the vector field as the electromagnetic field, 
and we discuss this model only as a mechanism for producing non-Gaussianity.^ Indeed the 
vector modes are coupled to the inflaton field by the same {(j)) F'^ term that generates 
them, and they source inflaton perturbations through this coupling. These perturbations 
add up incoherently to the inflaton vacuum modes, and are highly non-Gaussian. 

Ref. [78] computed the resulting bispectrum, (CkiCk2Ck3), in the equilateral configura- 
tion, i.e., ki = k2 = k^. The full bispectrum was computed in Refs. [24, 25]. The computa- 
tions of Refs. [24, 25] are restricted to n = 2 and n = —2 which produce, respectively, scale 
invariant "magnetic" and "electric" perturbations. The model enjoys a symmetry / -H- j, or 
n •H- — n, under which \E\ -H- Both zbn result in the same equation for ^. For brevity of 
exposition, we only refer to the n = —2 case in the reminder of this subsection. In this case 
one obtains, at the leading order in slow-roll, 

for the mode functions of each of the two polarizations of the "electric" and the "magnetic" 
fields in the super- horizon regime [25]. We note that the power in the "electric" field is frozen 
outside the horizon and scale invariant, whereas the power in the "magnetic" field decreases 
to negligible values. 

Let us assume that, at the beginning of inflation, say at the time t = 0, the "electric" field 
has a classical homogeneous value, E'', all across the universe with negligible perturbations. 
For n = —2, the classical equations of motion for the vector field are solved by a constant, 
E = E*^. This quantity is, however, not the classical "electric" field that would be measured 
by a local observer at t > 0, which we denote by E^. In fact, the modes given by Eq. (2.10) 
become classical after they leave the horizon (we denote them as infra-red (IR) modes), and 
they add up with E'' to give Ed. A given IR mode of wavelength A averages to zero on 

'""These problems persist also for a general evolution of I beyond the a" scaling [46]. 

^We continue to adopt the "electromagnetic" decomposition (2.8) for notational convenience. 
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regions of size L ^ X, but it is constant in each region of the size L <C A, and it adds up 
stochastically with E*^ and with all the other modes with A S> L generated during inflation 
to determine the value Ed in that region. An observer at time t > during inflation can 
only experience the value of Ed in its local Hubble patch. The average measured by this 
observer is drawn from a Gaussian distribution with the mean E'' and the variance given by 

E-E = 5-/ —k^El = — ttN, (2.11 

^ ^ {2'nf JHa{t=o) k 2n^ ' ^ ' 

where N is the number of e-folds from the start of inflation to the time t. The lower limit 
in the integral corresponds to the modes that left the horizon at the start of inflation (larger 
modes are not generated), while the upper limit corresponds to the modes that left the 
horizon at the time t (larger-momentum modes are still in the quantum regime and do not 
contribute to the classical average, Ed, in the Hubble patch of length ■^jjjjj)- 

The situation is completely identical to what happens to the so-called "stochastic infla- 
tion [79]," in which the variance of a massless scalar field, %, is determined by the stochastic 
addition of the IR modes, and grows as (x^) oc H'^N during inflation. It is well established in 
that context that this variance contributes to the theoretical expectation value of the scalar 
field measure by local observers. This is customary used, for instance, in the Affleck-Dine 
model of baryogenesis [80] or in the curvaton field [81]. The fact that the vector field has spin 
1 does not make any difference for these considerations, which simply follow from Eq. (2.10). 

Let us consider a mode, Cfej of a given comoving momentum, k. This mode leaves the 
horizon N/^ e-folds before the end of inflation. We are interested in the modes that affect the 
CMB anisotropies. Such modes leave the horizon A^^^. ~ A'cmb — 60 e-folds before the end of 
inflation. The Hubble patch that they exit is the one that eventually becomes our Hubble 
patch. When the mode Cfc leaves the horizon, the classical average of the "electric" field, 
Ed, in this Hubble patch is drawn from a Gaussian distribution with the mean E° and the 
variance given by Eq. (2.11), with A^ = Atot — A^fc, where A^tot is the total number of e-folds 
of inflation. 

In the presence of this mean, the kinetic term given in Eq. (2.7) results in the coupling 
of £int — 4a^Ed • 5E(^. This is the dominant operator for the part of C sourced by the vector 
field [25]. The power spectrum of C generated by this mechanism is given by 



Pc(fc) = Pf^(fe)[l + 5*(A:)cos^%,^J , gAk)^-^^f^^ , (2.12) 

where P^^^ = ^ sttCm^ square amplitude of the standard vacuum modes. The second 

term in Eq. (2.12) is the contribution of the sourced part of which (as phenomenologically 
required) we have assumed to be subdominant.^ It follows from the A^| proportionality 
that this term continues to grow in the super-horizon regime. This power spectrum was 
previously obtained in refs. [82-84] in the context of the anisotropic infiationary model [85], 
where however ^'ci was identified with E'^, missing the IR contribution. 



^Non-detection of statistical anisotropy in the WMAP 9-year data after the correction of non-circular beam 
effects [13] would imply a conservative upper bound of g, < 0.1. 
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The 5E mixing results in the following bispectrum in the squeezed limit, k-^ <^ ki ^ k2 



B^{ki,k2M) - '^^P^^^ iki)Pf^ (ks) \g* (ki) \Nk, 

The predicted power spectrum (Eq. 2.12) and bispectrum (Eq. 2.13) break statistical 
isotropy, as Ed picks out a preferred direction. However, the prediction for an isotropic 
measurement is obtained by averaging Eq. (2.13) over all directions of Eci- We find 

B^{ki,k2, ks) ~ 8P.(°) (ki) ) (ks) \g, (ki) \Nk, (1 + ^23) , (2.14) 

isotropic measurement ^ ^ 

where /uia = ki • ka. From this we obtain the Legendre coefficients in Eq. (1.1) as 

0.1 60 2 ^ ^ 

Due to simplicity of the model, Eq. (2.15) is a very predictive result, relating the bis- 
pectrum coefficients to the amount of statistical anisotropy of the power spectrum, i.e., 5*. 
This result holds for all models characterized by Eq. (2.7) and scale invariant "magnetic" 
or "electric" modes, including many analyses of the magnetogenesis mechanism [72] and of 
anisotropic inflation [85] (for which the departure from scale invariance is negligibly small; 
we note that in this work E*^ evolves on an attractor solution, but this has no consequence 
for the accumulation of the IR modes). An analogous result will also hold for the model of 
Ref. [86] and for the mechanism of Ref. [87], for which the scalar perturbations have been 
studied in Ref. [88]. 

The smallness of (7* limits the level of non-Gaussianity. However, a larger bispectrum, 
for a given value of 5*, can be obtained if the model is more complicated. For instance, 
one can arrange for a triplet of U(l) vectors, and assume that they have classical vacuum 
expectation values which are orthogonal to one another and of equal magnitudes [89]. In 
this case the power spectrum is statistically isotropic (5* = 0). This requires to assume that 
the IR sum is subdominant, as there is no reason to assume that the IR modes of the three 
vectors add up to orthonormal values. A larger bispectrum can also be obtained if there 
are additional fields and additional couplings, as in the waterfall mechanism of Ref. [90], in 
which a vector field of the kinetic term given in Eq. (2.7) is also coupled to the field that 
determines the end of hybrid inflation (see Ref. [25] for more detailed discussion) . 

2.3 Solid inflation 

Ref. [26] studied a rather unusual model, in which inflation is driven by a system which has 
a field-theoretical description of a solid. An equivalent version of the model was proposed by 
Ref. [91], under the name of "elastic inflation." 

Each volume element of the solid is characterized by a comoving label, (for instance, 
it can be the position of that element at the initial time t = 0). The functions, (i,x). 



*See Ref. [25] for the full expression; due to different conventions, the bispectrum of Ref. [25] is the one 
given here divided by {2t:)^^^ . 
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specify which volume element is located at a given position, x, at a given time, t. A solid at 
rest in comoving coordinates then obeys 



{<P') = X , = y , (03) = ^ (2.16) 

or, in short, {cj)^) = x*. Even if the vacuum expectation value of each field is x-dependent, 
a homogeneous and isotropic Friedmann-Robertson- Walker solution can still be obtained by 
requiring that the Lagrangian that controls the solid be invariant under rigid translations, 
0* + a\ and S0(3) rotations, (p'' — ^ Ojt/)^. 

At the lowest order in a derivative expansion, the translational invariance is guaranteed 
by considering functions 

B'^ =g^^'^d^<l,^d,<l>i, (2.17) 

and isotropy is obtained by requiring that the Lagrangian is a function of S0(3) invariants 
built from B'^K Only three independent invariants exist, and Ref. [26] chose 



5soiid = / d^x^F [X,Y,Z] , X = [B] , y = , Z = i-^ , (2.18) 



[Bf ' [BY 



where the square parenthesis denotes the trace of the corresponding matrix, e.g., [B] = 
The system has the energy-momentum tensor, T'^ = dia,g{—p,p,p,p), with [26] 

p = -F , p = F-\fx , (2.19) 

where the subscript denotes a partial derivative, and F and Fx are evaluated on the back- 
ground solutions given hy X = Y = and Z = ^. These invariants are chosen in 
such a way that X is the only one affected by the overall physical volume expansion; this 
immediately explains as to why only the derivative of F with respect to X enters into the 
expression for the pressure. 

Inflation is possible only if F is only mildly affected by the physical expansion, or, 
equivalently, only if Fx is sufficiently small. Specifically, we require 

e^-§ = ^«l. (2.20) 

One also requires Fxx to be small, so that r/ = ^ <C 1 [26]. 

Let us now discuss cosmological perturbations in this system. It is convenient to work 
in a spatially flat gauge, where the dynamical scalar and vector perturbations are all encoded 
in the perturbations of the scalar fields: 

S(l)' = Tt' {t, X) = TTL + TT^r, (2.21) 

where the vector components are transverse, diir^ = 0. The scalar and vector modes are 
decoupled from each other at the linearized level. ^ At the lowest order in the slow roll 



®As always in cosmology, the sectors of scalar and vector perturbations also include non-dynamical modes 
which are, in the spatially flat gauge, encoded in the Sgofi metric components. These fields can be integrated 
out as explained in [26]. This affects the action for the dynamical modes, ttl and tt^, at long wavelengths. 
Finally, there are also tensor perturbations - the gravity waves - encoded in the Sgij metric components, and 
which we do not discuss here. 
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parameters, and in the deep sub-horizon regime, the sound speeds of the scalar (longitudinal) 
modes, cl, and vector (transverse) modes, ct, are given, respectively, by [26] 

SO that the propagation is subluminal and non-tachyonic (0 < < 1 and < < 1) for < 

Fy + Fz < |e|F| [26]. Namely, the requirement of an accelerated expansion forces Fx to be 

small, while subluminality also requires that the combination Fy + Fz be small. Finally, the 

theory involves derivative interactions of the "phonon" fields, vr*, which necessarily become 

strong at some scale A. A detailed study in Ref. [26] gives an estimate of A ^ H (so 

that the linearized theory is also valid in the sub-horizon regime, up to A), provided that 
„ / „ N 2/3 

ec^ ^> ijj^j 1 which can always be obtained for a sufficiently small H. 

In a conformally flat gauge, the gauge-invariant scalar perturbation, C, evaluates to 
= —H^ = ^dir. Its solution exhibits two features that are peculiar in scalar-field infla- 
tion models, but that were nevertheless already seen in the models studied in the previous 
subsection [24]. The first one is the fact that ^ is not conserved on super-horizon scales, due 
to the anisotropic stress that does not vanish on super-horizon scales [26]. Indeed, following 
Ref. [26], we obtain 

STij,sc^Ur = a^M^HC [2 (3 - 2e + ??) - (3 + 3ci - 2e + r/) (sh % - J^j)] . (2.23) 

Let us recall that, also for the model described by Eq. (2.7), the anisotropic component 
of the stress-energy tensor, oc Ed^iSEj, does not vanish outside the horizon, and thus the 
anisotropic term in Eq. (2.12) grows outside the horizon. 

The second feature is that, analogously to the model described by Eq. (2.7), the bis- 
pectrum of solid inflation is largest in the squeezed conflgurations, and exhibits a nontrivial 
dependence on the angle between the modes in the squeezed limit. The dominant contribu- 
tion to the bispectrum is given by the interactions of vr encoded in the scalar field Lagrangian, 
while the metric perturbations provide a negligible contribution. The dominant interaction 
in a slow roll expansion is given by [26] 



r 



81 9 9 3 



(2.24) 



Also in this respect, the situation is analogous to the model discussed in the previous sub- 
section, where the dominant contribution to the bispectrum is obtained from Eq. (2.7), 
disregarding metric perturbations. When written in terms of this interaction exhibits non- 
trivial dependence on the direction of the modes which does not vanish in the squeezed limit, 
imprinting the nontrivial angular dependence in the bispectrum. 

At the leading order, the Legendre coefficients in Eq. (1.1) are given by [26] 



co~0 , C2 = 0(l)^ — . (2.25) 



Fy 1 

-L 

Namely, in the squeezed limit, the dominant contribution is given by the quadrupole term, 
whereas the monopole term is negligible. The dominant quadrupole term is essentially pro- 
portional to a free combination of parameters (we recall that avoiding superluminality and 
strong coupling at p < H imposes restrictions on the combination Fy + Fz but not on Fy 
or Fz individually). 
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3 Signatures in the cosmic microwave background 



In this section, we shall derive the flat-sky (Sec 3.1) and full-sky (Sec 3.2) formulae for the 
bispectrum of CMB temperature anisotropy from the bispectrum of curvature perturbations 
given in Eq. (1.1). We then calculate, in Sec. 3.3, the error bars of cq, ci, and C2 expected 
for a cosmic-variance-limited experiment measuring temperature anisotropy up to £ = 2000. 



3.1 Flat-sky formula 

While the full-sky formula is eventually needed for the analysis of full-sky temperature maps, 
let us derive first the flat-sky formula, as the flat-sky formula is usually simpler and more 
intuitively understandable. 

Under the flat-sky approximation, which is valid only on sufficiently small angular scales, 
i ^ 1, CMB fluctuations on the sky are expanded using the two-dimensional Fourier trans- 
form, instead of the spherical harmonics. The Fourier coefficients of CMB anisotropy, a{i) 
are related to the three-dimensional coefficients of the curvature perturbation, C(k), as [92] 



a{e) = £' dr 1^ (kll = Sj [k = ^kl + {l/D)\r) -^e^^^^ , (3.1) 

where k" = (fe^, ky), D = tq — t is the conformal distance out to a given epoch r, tq is the 
present-day conformal time, and Si is the so-called source function. ^'^ This relation simply 
tells us that a{£) measures the (" modes that are perpendicular to the line-of-sight direction 
(i.e., the modes on the sky), and the line-of-sight modes are washed out by integration. 
It is straightforward to compute the bispectrum of a(£) following, e.g., Ref. [92]: 



{a{ei)a{£2)a{£3)) = (2^)2<^(2) ^ ^ ^ cl6^(£i, £2, ^3) , (3.2) 



where 



6^(^l,^2,^3) 



r^dr 



n 

.n=l 



dT„ 



'2^ 



.n=0 



'1 • 



P^iki)P^ik2) + (2 perm) , 



(3.3) 



with 



g{e,k,T,r) 



'12 

kD 



-1/2 



Sl{k,T)j^ COS 



^"The source function is related to the radiation transfer function defined in Eq. (3.9) as Te{k) — 
dr Siik,r)j,ikD). 



The other kernel functions, J'^\ for L < 2 are given by 



r(0) 



1, 



(1) 



n 



Ft 



,(2) 



n 



n=l V fc„D. 



and 



•'I 



(0) 



Fl 



(0) 



(1) 




kn{r - Dn) 



n=l k„D„ ' 

(3.5) 
(3.6) 
(3.7) 



While these formulae are still complicated, one can read off the leading behaviours 
of these expressions in the small-scale limit, in which the dominant contributions in the k 
integration come from the modes with k ~ i/D. In this limit the kernel functions become 



•^1 



(0) 



0, Fr 



(1) 



0,FI 



(1) 



l,F 



(0) 



1/2, and ^ 3/2. We thus find 





3 „ 


/ r^dr 


n / 


1 — oo 


n=l 



xPi(£i • i2)Pdh)Pd>^2) + (2 perm). 



dk 

~^ Q{^ni kni Tji, r) 
ZTT 



(3.8) 



This result shows that the CMB bispectrum is proportional to Pii^i ■ ^2) (and its permuta- 
tions), which is expected from PlO^i ■ ^^2) (and its permutations) in the three-dimensional 
bispectrum of the primordial curvature perturbation. 

While the approximation of ^ ~ kD gives a simple and transparent result given in 
Eq. (3.8), it is less precise than the original form given by Eq. (3.3). In Appendix A, we 
discuss the precision of Eqs. (3.3) and (3.8) with respect to the full-sky result given in the 
next subsection. 

3.2 Full-sky formula 

Encouraged by the flat-sky results, we now move onto the full-sky case. The CMB tem- 
perature anisotropy on the celestial sphere is expanded by the spherical harmonic function 
as (5T(n)/r = m "^m^^m(n), where fi is a three-dimensional unit vector pointing toward 
a given direction on the sky. The spherical harmonics coefficients, aim, are related to the 
primordial curvature perturbation as 

d^k 



4Tr{-iY 



(2^)3 
edk 



r,(A:)Ckl^;„(k) 
r^Te{k)Cim{k), 



(3.9) 



where 7e{k) is the so-called radiation transfer function, and we have defined the curvature 
perturbation expanded in spherical harmonics as Ctm{k) = / (i^k (^k^m(^)- Then, the bis- 
pectrum of aim can be straightforwardly calculated as 




n47r( 



.n=l 



k'^dkn 
"(2^ 



Te„ {k„ 



n C£nm„(fc„ 



(3.10) 



\?i=i 
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with the bispectrum of dm related to B(^{ki, ^2,^3) as 



\n=l 



3 

.n=l 



(27r)35(3) (ki + k2 + kg) ^2, A:3).(3.11) 



Using the bispectrum of C given in Eq. (1.1), B(^{ki, k2, ks) can be expanded as 



47r 



B^{ki,k2,k3) = Pdh)Pdk2)^CL^^^Yl,j{ki)YLM{k) + (2 perm). (3.12) 

L M 

Using the definition of the delta function, (5^^^(k) = (27r)~^ J d^x e^^'^, we also expand the 
delta function as 



5(3) (ki + k2 + kg) = 8 / 

Jo 



r'^dr 



n E (-l)'^i^n(A:nr)y£„M„(kn) 



.n=l LnMn 

Li L2 



lLiL2L:j , 



Ml M2 M3 

where the 2x3 matrix denotes the Wigner-3j symbol, and the / symbol is defined by 



(3.13) 



(2/i + l)(2/2 + l)(2/3 + l) fh h h 





47r 



(3.14) 



Now, performing the integrals of the spherical harmonics over ki,k2, and k3, and per- 
forming the summations over Mi, M2, M3, and M as described in Ref. [93], we obtain 



{2TTfB(^^i,^i^jdki,k2,kz) 



\n=l 



ti l2 ^3 

mi 1712 



(3.15) 



where 



-Sc/i^2^3(^i'^2, A;3) = 8 / r'^dr 

Jo 



.n.=l Ln 

4tt 



IL1L2L3 



2L + 1 



h ^2 k 
L2 Li L 



Pdki)Pdk2) + (2 perm) 



(3.16) 



Here, the 2x3 matrix enclosed by curly brackets denotes the Wigner-6j symbol. As the 
primordial bispectrum given by Eq. (1.1) is rotationally invariant, the bispectrum expanded 
in spherical harmonics must also be rotationally invariant. This means that the dependence 
of the bispectrum on mi,m2 and must be given by the Wigner-3j symbol, as shown in 
Eq. (3.15). This property ensures rotational invariance of the CMB bispectrum. 

Substituting Eqs. (3.15) and (3.16) into Eq. (3.10), we finally obtain the full-sky formula 
for the CMB bispectrum: 



n 

\n=l 



Hi l2 *3 

mi m2 



h i2 h 
mi 7712 7713 



(3.17) 
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where 

POO ^ 



.n=l L„ 



/£iL,l/..l.l(-1)'^+^Ml3/3 { l\ l\ } + (2 perm) , (3.18) 



2L + 1 
and 



a,(r) = - / k^dk %{k)n{kr) , (3.19) 
71" Jo 

/3,i(r) = - / fc^^A: Pcik)Ti{k)jL{kr) . (3.20) 



Owing to the selection rules of the Wigner symbols, ii, (.2 and £3 are constrained by parity 
invariance and the triangular condition: 

^1+^2 + ^3 = even , \h - l2\ < h < h + ^2 ■ (3.21) 

The former constraint is a consequence of the bispectrum of curvature perturbations given 
by Eq. (1.1) being parity-even. The summation ranges of Li and L2 are also restricted to 

Ln = \in-L\Mn-L\+2,--- , £n + L - 2, 1^ + L . (3.22) 

In the full-sky formula given by Eq. (3.18), the angle dependence for L > induces 
a coupling among £i, £2 and £3 via the Wigner-6j symbol. As a result, Eq. (3.18) is not 
separable (or at least not obviously separable) with respect to ts unlike the usual local-form 
CMB bispectrum without angle dependence, i.e., L = 0. 

Figures 2 and 3 show the absolute values of the full-sky reduced bispectra, bf^i^^^ = 

^tii2i3 ^^f-i^ih)"^ ^ = 0) 1) ™d 2. Note that the full-sky reduced bispectrum reduces 

to the flat-sky bispectrum, b{ii,i2,£3), that we discussed in the previous subsection, in the 
small-sky limit [7]. 

Figure 2 shows the equilateral triangles with i = ii = £2 = while Figure 3 shows 
triangles with £1 = £2 = 200, which become squeezed triangles for £3 <C 200. We find that 
the amplitudes of the equilateral triangles monotonically decrease as L increases. We can 
understand this by using the flat-sky formula given in Eq. (3.8): the Legendre polynomials 
give the ratio of L = 0, 1, and 2 terms as 

b'^=\£,£,£) : b'^=H£,£,£) : b'^=\£,£,£) = 1 : -1 : -i, (3.23) 

for £,.£. = -I (^i^j). 

Figure 3 shows the squeezed triangles with £3 <^ £1 = £2 = 200. In the squeezed limit, 
the CMB bispectrum of L = 1 is highly suppressed compared with those of L = and 2. 
This is simply due to symmetry: the L = 1 term vanishes in the exact squeezed limit. Again, 
the flat-sky formula given in Eq. (3.8) gives the ratio of L = 0, 1, and 2 terms as 

6^=0(£i,£i,4) : 5^=H£i,4,^3) : 6^=2(£i, 4, ^3) = 1:0:-^, (3.24) 

for ii-is = = £2 •4- 

As these calculations are quite involved, we provide the simplified analytical formula 
in the Sachs- Wolfe limit in Appendix B. This test validates our numerical results shown in 
Figures 2 and 3. 



- 14 - 



0.1 F 



X 



^ 0.0001 r 



+ 



le-05 r 



le-06 r 



le-07 




L = 
L = l 
L = 2 



_i I I ' ' ■ ' 



10 



100 



1000 



Figure 2. Absolute values of the equilateral CMB temperature reduced bispectra, for L = Q 

(solid), L = 1 (long-dashed), and L ~ 2 (short-dashed). 



3.3 Expected uncertainties on ci and C2 

In this subsection, we calculate the l-a error bars of cq, ci and C2, i.e., 5co, 5ci, and 5c2, 
expected for a cosmic-variance-limited experiment measuring temperature anisotropy. Here, 
we shall focus on a simultaneous estimation of a pair of parameters: (cq, ci) and (cq, C2). We 
give the full constraint varying all three parameters simultaneously in Appendix C. 
Following Ref. [7], we calculate the Fisher matrix, -Fll', from 



2<£l<<?2<<;3<£max 



a 



(3.25) 



where the variance of the CMB bispectrum, cf^^^^s ' given by 



W2^3 - Cl^Cl^Cl^ (-1) 



(3.26) 



with Ci being the power spectrum of temperature fluctuations. As we consider a cosmic- 
variance-limited experiment, we ignore instrumental noise here. 

As we show in Appendix C, cq and ci are nearly uncorrelated, so are cq and C2; however, 
ci and C2 are highly correlated. Therefore, in this subsection, we shall consider submatrices 
of -Fll' involving only either (co,ci) or (co,C2), and study the full matrix in Appendix C. 
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Figure 3. Same as Figure 3, but for tlie squeezed triangles, \bi^i^i^\, witli li ~ t2 = 200, as a function 



We define the submatrix (a 2 x 2 matrix) as 



(2)_p. 



-foo Fql 
Flo Fll 



(3.27) 



wliere L takes on either 1 or 2. The 1-a marginalized error bars are then given by the matrix 
inverse as {5co,6cl) = 



(2) Z?-l 



In Figure 4, we show the ratios of error bars, 6ci/5co and 602/ 6co, as a function of 
the maximum multipole in the sum, ^max- The solid and short-dashed lines show the exact 
results for L = 1 and L = 2, respectively, while the long-dashed and dotted lines show the 
corresponding Sachs-Wolfe approximations for L = 1 and L = 2, respectively. We find that 
the Sachs- Wolfe approximations trace the overall behavior of the exact calculations well. 

The error bar on ci is an order of magnitude larger than that on cq for imax ^ 100, as 
the L = 1 bispectrum has a vanishing amplitude in the squeezed limit. On the other hand, 
the error bar on C2 is comparable to that on cq: the Sachs- Wolfe approximation gives an 
asymptotic relation of 5c2 = 4:6cq. The exact calculation gives 6c2 ~ SScq for irna.^ = 2000. 

Finally, the 1-a error bars expected for a cosmic-variance-limited experiment measuring 
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Figure 4. Ratios of tlic expected error bars, Sch/Sct) (L = 1 and 2), as a function of the maximum 
multipoles in the sum, ^max- The solid and short-dashed lines show the exact results for L — 1 
and L = 2, respectively, while the long-dashed and dotted lines show the corresponding Sachs- Wolfe 
approximations for L = 1 and L = 2, respectively. We find that the Sachs- Wolfe approximations 
trace the overall behavior of the exact calculations well. 



temperature anisotropy up to £max = 2000 are given by 



(<5co,5ci) = (4.4,61) , (3.28) 
(5co,(5c2) = (4.4,13) . (3.29) 



See Eq. (C.l) for the full Fisher matrix. 



4 Consistency relations with higher spin fields 

Primordial correlation functions in the limit that some combinations of external momenta go 
to zero - soft limits - play a special role in constraining the physics of inflation. Significant 
squeezed non-Gaussianity is associated with the presence of extra light degrees of freedom 
during inflation, hence soft limits can be understood as probing the spectrum of light fields 
in the early universe. Moreover, soft limits are observationally relevant and are subject to a 
number of interesting theoretical consistency relations. The first example of such a consis- 
tency relation was noted in [12] and established under much more general conditions in [1]: 
limfc3^o ^c(^i; ^2; ^3) = (1 ~ 'T's)-Pf (^i)-Pf (^3)- This holds independently of the inflationary 
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Lagrangian, under the assumption that there is only a single field, an attractor solution has 
been reached [5, 6], and the initial state is in a Bunch-Davies state [3, 4]. 

Our new parametrization given by Eq. (1.1) represents a non-trivial modification of this 
consistency relation: 



lim Bc{ki,k2,k3) 



clPl{^i ■ ks)^ Pc(^i)^c(^3) . (4.1) 



The possibility that some of the cl coefficients can be ^ |n<j — 1| ~ 10~^ indicates extra 
light fields in the early universe, while the non-trivial angular dependence is associated with 
anisotropic sources (such as higher spin fields). We also expect non-trivial soft limits for 
higher-order correlation functions. To explore such effects, we introduce the quantities 



^.^ 5 B(^{ki,k2,k'i) 
fc3™o 12 P^{ki)P^{k3) ' 

1 T(;{ki,k2,k3,k4) 



lim 



ki+k2^o 4 (|ki + k2|)Pc(A:i)i^c(^3) 



(4.2) 
(4.3) 



where (CkiCk2Ck3Ck4) = (2vr)^5(^)(^j kj)T^(A;i, ^2, ^3, ^4). For the local-type non-Gaussianity 
the quantities and become the standard (momentum independent) non-linearity 
parameters. In this case there is an interesting consistency relation: 

2 



TNL ^ [ T^JNL 



(4.4) 



where we dropped the superscript "eff" to emphasize that this relation is understood in the 
case where Eqs. (4.2) and (4.3) are independent of momenta. This relation was first noted 
by Suyama and Yamaguchi in [16] and further explored in [17-21]. 

In this subsection, we explore the non-Gaussianity consistency relations in the context 
of the /2(^)f2 model given in Eq. (2.7). In general this model breaks statistical isotropy, as 
discussed in Sec. 2.2. From Eq. (2.13) we obtain 



-feff ^ 1 n Ar2 AT ^'^-^cl 
/^^^107V,^iV,3-^ 



cos^e 



ki , E^i 



COS^ 01. p + COS Ol 

ks,Eci ki, 



p COS Of. A COS 9 1 t 



(4.5) 

where we recall that \g*{Nk^)\ = N^, tV{4>) ^ assumed. 

We compute the trispectrum for the first time in the model given by Eq. (2.7). The 
computation follows very closely the analogous one of the bispectrum performed in [25], and 
so we omit the technical details. We obtain 



^cff ^ 



144iV|.iV|3 '''''' 



■■eVi<P) 



1 - cos^ 



ki, -Ed 



cos^ Or. 



+ cos 01 fp COS 9t A 



fcl2, -Eel 
COS " 



cos^ Oi 



k3,E^ 



fci, fca 



+ COS 9; fp COS 9t fp COS 9t i 

kl,Ecl Kl2,-Bcl Kl,fel2 

+ COS 9r A COS 9r cos 9r ,- 

Ki2,-Eci k3,Eci ki2,k3 



COS 9; £, COS 9f A COS I 

ki,Ecl k3,Ecl 



h f COS^r f 



(4.6) 
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where ki2 is the unit vector in the direction of ki +k2. We note that the unit-vector ki2 enters 
in this expression, even though the corresponding vector ki + k2 vanishes in the squeezed 
hmit (analogously to the fcs-dependence of the bispectrum) . 

We observe that, in the l'^{(f>)F'^ model, both and exhibit highly non-trivial 
momentum dependence; thus, it is not sensible to compare different configurations. One can 
readily see that r^^ vanishes for several configurations. This happens, for example, if one of 
ki , ki2 , ^3 is parallel to Ed, while the other two vectors are perpendicular to Eci- We do not 
interpret this as a violation of the Suyama-Yamaguchi inequality; the expression given by 
Eq. (4.4) requires either that the non-linearity parameters are momentum-independent, or 
else that and (/^) have the same momentum dependence so that one can factor out 
an amplitude which obeys Eq. (4.4). The original form of the Suyama-Yamaguchi inequality 
is simply not applicable to the model given by Eq. (2.7).^^ 

In Ref. [21] Assassi et al derived a general inequality relating the soft limits of the 
bispectrum and trispectrum: 

j d^qiSq2 r^^i(qi,k-qi,q2,-q2-k) Pc(9i)^cfe) 



> 



/ A^/fL(q,-q-k,k)Pc('z) 



(4.7) 



where the A; — )• limit is understood. This inequality reduces to Eq. (4.4) in the local 
case, but is completely general and should be respected by any model. We can easily verify 
that, although the Suyama-Yamaguchi inequality is not meaningful here, Eq. (4.7) is still 
respected. After evaluating the angular integral, the left hand side of Eq. (4.7) becomes 

LHS^ 1024vr2-^sin2^^^^^^(^y rfgg2^,P^(g)j , (4.8) 



and, proceeding analogously, the right hand side becomes 
RHS ~ imAii'^Nl 



^cl 



Therefore, we find 

LHS 



RHS (iVfc) |sin^% ^ 



sin^^ ^^^^d dqq'NgP^iq)^ . (4.9) 



> 1. (4.10) 



Note that <C 1 has been assumed (as also required by phenomenology) throughout this 
subsection. We therefore conclude that the integrated inequality given by Eq. (4.7) is satisfied 
for any orientations of k relative to the classical vector field background. 



^^Ref. [94] presents a detailed discussion on how to generalize the Suyama-Yamaguchi inequality to 
momentum-dependent non-linearity parameters, with a particular emphasis on the case of broken statisti- 
cal isotropy. The fact that we have shown that there exist some configurations with vanishing implies 
that, at least in principle, one may obtain observational evidence for a smaller squeezed trispectrum than the 
one obtained from scalar fields, provided that one can construct an observable quantity which is sensitive to 
those configurations. It remains to be seen whether such a measurement is feasible. 
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5 Conclusion 



The angle dependence of the bispectrum of primordial curvature perturbations in the squeezed 
configuration is sensitive to the presence of vector fields and non-trivial symmetry dur- 
ing inflation. In this paper, we have explored phenomenological consequences of the new 
parametrization of the bispectrum given by Eq. (1.1): i?( = "^j^ CiPL{ki ■ k.2)P({ki)P(^{k2) + 
(2 perm). This form is physically well motivated, and we have given three examples in Sec. 2: 
the curvature perturbation sourced by the anisotropic stress of magnetic fields; that sourced 
by an interaction with a vector field of the form P((j))F'^; and solid inflation. 

We find that a cosmic- variance- limited CMB experiment measuring temperature anisotropy 
up to ^max = 2000, such as the Planck satellite, can measure ci and C2 down to 6ci = 61 and 
6c2 = 13 (68% CL). The latter error bar is comparable to (and only a factor of three larger 
than) the error bar of cq = 6/nl/5; thus, if the forthcoming Planck data reveal evidence 
for Co, one should also measure C2 to understand the nature of sources of non-Gaussianity. 
Moreover, even if the Planck data do not reveal evidence for cq, one should still measure C2, 
as solid inflation can generate large C2 without generating detectable cq. Sensitivity to ci is an 
order of magnitude worse than that to cq or C2 because the term proportional to ci vanishes 
in the squeezed limit due to symmetry. 

The angle-dependent bispectrum in the squeezed configuration is a natural consequence 
of broken statistical isotropy. Broken isotropy also leads to a non-trivial modification of 
the inequality between the local-form trispectrum amplitude, tnlj and /^l- We find that 
the original form of the Suyama-Yamaguchi inequality, tnl > (6/nl/5)^, does not apply 
to the current model, due to the momentum- and shape-dependence of tnl and /nl- For 
example, we find some squeezed configurations in which tnl vanishes. It remains to be seen 
how sensitive the forthcoming tests of the Suyama-Yamaguchi inequality using the Planck 
or the large-scale structure data are to the decrease of the trispectrum amplitudes for these 
particular shapes. We also find that a general inequality of Ref. [21] is satisfied in this model. 

What is next? Phenomenological consequences of Eq. (1.1) for large-scale structure 
of the universe such as the dark matter halo bias [95], bispectrum, and trispectrum should 
certainly be explored. For instance, a consistency relation [96, 97] between the squeezed 
bispectrum and the power spectrum of dark matter density fluctuations has been proved 
at the full nonlinear level, and with an initial isotropic non-Gaussianity, in [96]. It may be 
interesting to study a signature of anisotropic initial non-Gaussianity in that context. Also, it 
is quite possible that the coefficients cl depend on wavenumbers, which would be particularly 
interesting for dipolar and quadrupolar modulations of the observed power spectrum in our 
sky [33]. Indeed, in the case of the /^(0)F^ model discussed in Sec. 2.2, the cl coefficients 
exhibit a logarithmic running with wavenumber. It would be interesting to study possible 
effects of such wavenumber dependence. 
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Figure 5. Absolute values of the CMB temperature redueed bispeetra. The solid, long-dashed and 
short-dashed lines show the full-sky results for L = 0, 1, and 2, respectively, while the plus, cross, and 
star symbols show the simplified fiat-sky results from Eq. (3.8) for i = 0, 1, and 2, respectively. The 
square symbols show the original form of the flat-sky result for L = 2 from Eq. (3.3) before further 
approximation. (Left panel) Equilateral triangles, (Right panel) Squeezed triangles, 

with £i = £2 ~ 200; as a function of ^3. 
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A Precision of the fiat-sky approximation 

How precise are the flat-sky formulae given by Eqs. (3.3) and (3.8)? In Figure 5, we compare 
the full-sky results with the flat-sky results. We find that, for L = and 1, the simplified 
flat-sky formula given by Eq. (3.8) yields the bispeetra in the equilateral and squeezed config- 
urations which are in good agreement with the full-sky results at £ > 100. However, we flnd 
that, for L = 2, the simplifled formula systematically underestimate the magnitude of the 
bispeetra in both conflgurations. The equilateral result suggests that the simplifled formula 
provides an adequate result only at ^ > 800. 

While these results appear to suggest that the precision of the simplifled formula de- 
grades as L increases, this is not the case: the flat-sky results in the Sachs- Wolfe limit (which 
are not shown in this paper) show that the simplified formula overestimates the magnitudes 
of the bispeetra of L = 1 and 2 in the Sachs- Wolfe limit by a similar amount in both equilat- 
eral and squeezed configurations. Therefore, we conclude that the simplified formula given 
by Eq. (3.8) should only be used for quantitative calculations of L = or for qualitative 
calculations of L = 1 and 2, and the original formula given by Eq. (3.3) should be used for 
quantitative calculations of L = 1 and 2. Needless to say, the full-sky formula should always 
be used for the calculations involving multipoles of £ < 100. 
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Figure 6. Absolute values of the CMB temperature reduced bispectra for L = (solid), 1 (long- 
dashed), and 2 (short-dashed), in the Sachs- Wolfe limit. (Left panel) Equilateral triangles, \biif\ 
(Right panel) Squeezed triangles, {bf^^^ij, with £i = ^2 = 200, as a function of £3. 



B Analysis in the Sachs- Wolfe limit 

As the calculations presented in Sec. 3.2 are quite involved, some appropriate approximations 
would be useful for understanding the analytical structures of the basic results. 

The Sachs-Wolfe limit, in which the radiation transfer function is given by 7i{k) — )• 
— ij£(fer*), provides such a convenient approximation. With this transfer function, a^(r) 
(Eq. 3.19) simplifies to a^(r) — )• —■^S{r — r^), where = tq — r* is the conformal distance 

to the last scattering surface. Similarly, for a scale- invariant spectrum of C, -Pc(^) 
/3i{r) (Eq. 3.20) becomes 



2n^ 



A.. 



10 ^r(^^)r 



(B.l) 



where r(x) is the Gamma function. Using these and in Eq. (3.18), one finds the 
Sachs-Wolfe approximation of the CMB bispectrum as 
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Figure 6 shows the reduced CMB temperature bispectra in the Sachs-Wolfe limit. The 
basic behaviors, such as the monotonic decrease of the equilateral amplitudes as a function 
of L and the suppression of the L = \ term in the squeezed limit, are all reproduced by the 
simple Sachs-Wolfe limit calculations. These results may be compared with Figures 2 and 3. 

C Full Fisher matrix 



In Sec. 3.3, we have presented the l-fx marginalized constraints on cq, ci, and C2, assuming 
that only cq and one of ci and C2 are varied simultaneously. In this Appendix, we provide 
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the full Fisher matrix, F^l', involving all of cq, ci, and C2, calculated up to £max = 2000: 

/ 5232 16.94 -5.986 \ 
F^y = 16.94 26.53 66.85 x 10"^ . (C.l) 
V -5.986 66.85 618.1 / 

From this matrix, one can compute the cross-correlation coefficients, = / V FllFl'V ■ 
We find roi = 0.045 and ro2 = —0.003, indicating that cq and ci are nearly uncorrelated, so 
are cq and C2- However, there is a high degree of correlation between ci and C2'- = 0.522. 
As a result, the marginalized error bars increase slightly to 

(<^co, 5ci,5c2) = (Jf^\ \/f^, = (4.4, 72, 15) . (C.2) 
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